An Infrared Determination of the Reddening and Distance to 

Dwingeloo 1 

Valentin D. Ivanov, Almudena Alonso-Herrero, Marcia J. Rieke and Don McCarthy 
Steward Observatory, The University of Arizona, 933 N. Cherry Ave, Tucson, AZ 85721 

E-mail: vdivanov@as.arizona.edu, aalonso@as.arizona.edu, mrieke@as.arizona.edu, 

mccarthy@as.arizona.edu 



Received 



accepted 



- 2 - 



ABSTRACT 

We present for the first time infrared observations of the nearby highly 
obscured galaxy Dwingeloo 1 (Dwl), including deep if-band imaging covering 
a total of 4.9' x 4.9', together with J and K s imaging of the central 2.5' x 2.5'. 
We used the small dispersion of the intrinsic infrared colors of spiral galaxies to 
determine an infrared if -band extinction of Ah = 0.47 ± 0.11 mag towards Dwl. 
In using infrared colors, the uncertainties in the reddening and distance are 
reduced by a factor of three. The if-band magnitude corrected for extinction 
and the infrared Tully-Fisher relation are then used to estimate a distance 
modulus of (m - M) = 28.62 ± 0.27, and thus a distance of d = 5.3l^Mpc, 
which puts Dwl at the far end of the IC 342/Maffei 1 & 2 group. Our result is 
largely independent of the nature of the reddening law because we estimated 
both the reddening and the distance at the same wavelength range. 

Subject headings: galaxies: individual (Dwingeloo 1) — galaxies: distances and 
redshifts — galaxies: photometry — infrared: galaxies 
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1. Introduction 

Dwingeloo 1 (Dwl) is a large SBb/c galaxy, discovered both in a systematic Hi 
emission survey of the northern part of the Milky Way in search of obscured galaxies in 
the Zone of Avoidance by Kraan-Korteweg et al. (1994), and independently by Huchtmeier 
et al. (1995). The knowledge of the local mass distribution has implications for the 
peculiar velocity field, the direction and amplitude of the Local Group acceleration, the 
determination of parameters such as Qq and H , and on the understanding of the formation 
and evolution of groups of galaxies (e.g., Peebles 1994; Marinoni et al. 1998). The discovery 
of this galaxy proved a long standing suspicion that the tidal disruptions of Maffei 2 may 
be due to the presence of another massive galaxy nearby (e.g., Hurt et al. 1993). 

Dwl lies in the direction of the IC 342/Maffei 1 & 2 group of galaxies, about 2 degrees 
away from Maffei 2. This corresponds to a physical separation of 175 kpc assuming that 
Dwl and Maffei 2 are at a distance of 5 Mpc. Being the nearest barred spiral system, 
Dwl offers a unique possibility to study the effect of the bar at high spatial resolution. 
The discoverers classified Dwl as an SBb or SBc galaxy (T=4) and measured an angular 
diameter of 4.2arcmin. Later on, McCall & Buta (1997) re-classified it to SB(s)cd with an 
angular diameter of 9.9arcmin at \ii = 25.0 mag arcsec -2 based on deep optical /-band 
imaging. Burton et al. (1996) extensively studied the neutral hydrogen content of Dwl and 
measured Hi profile widths at 20% and 50% level of 201.2 ± 0.4km s^ 1 and 187.6 ± 0.6km 
s _1 respectively. The measured inclination of the gaseous disk was 51 ± 2 degrees and the 
position angle 112 degrees, with the major axis aligned with the bar. 

Since the discovery of Dwl the determination of its distance has been hampered by 
the poorly known Galactic extinction. Optical V RIH a imaging, long slit spectroscopy, and 
IRA S observations were summarized by Loan et al. (1996), who used a number of methods 
to estimate the foreground extinction towards Dwl. The optical color excesses yielded 
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A v = 7.8 ± 3.0 ,mag, the measured Galactic Hi column density A v = 4.5 mag, and the 
100 /im IRAS flux A v = 3.2 mag. Finally, they applied optical f and i?-band Tully-Fisher 
relations to obtain distances ranging from 1.3 to 6.7 Mpc, with an average value of about 
4Mpc (assuming H = 75km s -1 Mpc -1 ). Their main source of uncertainty was the value 
of the Galactic extinction. Phillipps & Davies (1997) challenged these extreme distance 
estimates on the basis of the very narrow span of central surface brightness in present 
day spiral galaxies. Their best estimate for the extinction (Ab — 6 mag) places Dwl at 
a distance of 3.1 — 3.6 Mpc. These authors also employed the diameter version of the 
Tully-Fisher relation (Persic, Salucci, & Stel 1996) and obtained a distance of 2.7 Mpc. 

The primary goal of this study was to put stronger constraints on the foreground 
extinction, and to obtain a better estimate of the distance to Dwl. We chose to use infrared 
colors because of their small intrinsic variations among spiral galaxies (Aaronson 1977). In 
addition, the extinction in the if -band is about three times smaller than in the /-band and 
about six times smaller than in the V-band (Rieke & Lebofsky 1985). Finally, the infrared 
Tully-Fisher (IRTF) relation shows a smaller intrinsic scatter than its optical counterpart 
(Aaronson, Huchra, & Mould 1979; Freedman 1990; Peletier & Willner 1993). The IRTF 
relation allows us to determine both the extinction and the distance at the same wavelength 
range, minimizing the errors arising from possible variations of the reddening law. 



2. Observations and Data reduction 

We obtained JHK S imaging of Dwl using a 256 x 256 NICMOS3 array at the 2.3-m 
Bok Telescope of the University of Arizona on Kitt Peak, with a plate scale of 0.6arcsec 
pixel -1 during a number of observing runs. We constructed a deep 4.9' x 4.9' i/-band 
mosaic of Dwl, whereas the J and K s images only covered the central 2.5' x 2.5'. Additional 
if -band imaging using a 1024 x 1024 array and plate scale 0.5 arcsec pixel -1 was obtained 
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at the same telescope on a subsequent observing run to calibrate the deep if -band imaging. 
The observational strategy consisted of taking galaxy images interleaved with sky images 
6' — 7' away from Dwl. Details of the observations are listed in Table 1. 

The data reduction included subtraction of dark current frames, flat-fielding with 
median combined empty sky frames, and sky subtraction. The mosaics were constructed by 
shifting the images to a common position with cubic spline interpolation. The photometric 
calibration was performed using observations of standard stars from the lists of Elias et 
al. (1982) and Hunt et al. (1998) when conditions were photometric. Conditions were 
non-photometric for the if-band mosaic, and this image was self-calibrated using the 
1024 x 1024 data. In the next section we will be making use of infrared colors of spiral 
galaxies to derive an estimate of the extinction to Dwl. The colors for spiral galaxies 
were obtained with the i^-band filter, whereas our measurements were taken with the K s 
(short-fT) filter. Therefore it is necessary to determine what the difference is between 
the two filters. Very recently Persson et al. (1998) have obtained a new set of JHKK S 
photometric standards. In their study the average difference between the K and K s 
magnitudes (for the red standards) is 0.0096 mag, with a standard deviation of 0.017 mag. 
We assume that the use of the K s filter instead of the K filter introduces an extra 
uncertainty in the photometric calibration of this filter of ±0.02 mag. The errors associated 
with the photometric calibration are 0.05, 0.07, 0.06 mag in J, H and K s respectively. 

In Figure 1 we display the iJ-band mosaic on a logarithmic scale. As can be seen from 
this figure, there is a large number of foreground stars which need to be removed prior 
to analyzing the data. The star removal from the JHK S images was done interactively. 
The affected pixels were then replaced with a linear surface fit to a circular annulus 
around each star. The automatic procedures failed largely because of PSF variations under 
non-photometric conditions. A bright star located southwest of the galaxy posed a major 
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problem, and was masked out throughout the data reduction. We assumed radial symmetry 
and replaced the region within about 90 arcsec from the star with the data from the 
opposite side of the galaxy. We performed the photometry on the cleaned images measuring 
the flux within elliptical isophotes with fixed position angle and ellipticity as determined 
from Hi observations (Burton et al. 1996). 



The surface brightness profiles in JHK S are given in Table 2, and displayed in Figure 2 
along with the radial distribution of the J — H , H — K s and J — K s colors, and the total 
apparent magnitudes in JHK S . In both Table 2 and Figure 2 the error bars represent the 
combined 3<r variations from photon statistics, sky background variations, the elliptical 
isophotal fitting, and the photometric calibration. From the radial distribution of the colors, 
it is clear that the center of the galaxy appears slightly redder than the outer regions, as it 
is the case with the infrared colors of most spiral galaxies (Terndrup et al. 1994; de Jong 
1996). Loan et al. (1996) report inverse optical color gradients in Dwl, however, this may 
be the result of active star formation along the bar as found in some barred spirals (Shaw 
et al. 1995). 

Before we tried to obtain an estimate of the extinction using the infrared colors, we 
fitted straight lines to assess the variation of the colors with increasing apertures, 



3. 



Discussion 



3.1. 



Colors and Extinction 



H-K s = 0.435(±0.042) - 0.072(±0.071)i? 



(1) 
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J-H = 1.099(±0.040) - 0.056(±0.067)i? 



(2) 



J-K s = 1.532(±0.042) - 0.125(±0.071)i? 



(3) 



where R is the semi-major axis in units of arcmin. These color gradients imply a 
change of 0.06 — 0.13 mag within the inner 2 arcmin of Dwl, which may cause significant 
uncertainties in the reddening estimate. The smaller field of view of the J and K s images 
prevents us from obtaining the total colors of Dwl. However, the total flux in the JHK S 
bands is dominated by the inner 2'-diameter region (see Figure 2, bottom panel). Hence, we 
adopt the total observed colors at radial distance of 1 arcmin from the center of the galaxy to 
be representative for the whole galaxy (J — H = 1.04 ± 0.10 mag, J — K s = 1.40 ±0.11 mag, 
H — K s = 0.36 ± 0.11 mag). To estimate the color excesses we compared these 
colors to the mean integrated colors of SBb-SBcd galaxies: J — H = 0.73 ± 0.02 mag, 
J — K s — 0.94 ± 0.03 mag, H - K s = 0.21 ± 0.02 mag (Aaronson 1977). We have increased 
the errors in Aaronson's colors to account for the uncertain Hubble type of Dwl. The color 
excesses were converted into if-band (A H ), visual (A v ) and 5-band (A B ) extinctions using 
Rieke & Lebofsky (1985, RL85) and Mathis (1990, M90) extinction laws. The results are 
summarized in Table 3. Henceforth we will use Ah = 0.47 ± 0.11 mag from Mathis (1990) 
extinction law for the sake of compatibility with previous work. This value is close to the 
estimate based on the IRAS 100 /im flux (Loan et al. 1996). 

Finally, we used the combined optical - near infrared colors to verify our result. Loan 
et al. (1996) reported the following total apparent (not corrected for reddening) magnitudes 
for Dw 1: m 7 = 10.7 ± 0.2 mag, m R = 12.2 ± 0.2 mag, and m v = 14.0 ± 0.5 mag. We 
measured a total if -band magnitude m H = 8.3 ± 0.2 mag, and compared the observed 
colors with the intrinsic colors as determined by de Jong (1996): I — H = 1.44 ± 0.20 mag, 
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R- H = 2.01 ± 0.20 mag, and V - H = 2.50 ± 0.20 mag. Mathis (1990) extinction law 
yields A H = 0.56 ± 0.23 mag, 0.58 ± 0.12 mag, and 0.47 ± 0.11 mag respectively, in good 
agreement with our infrared estimates. 

3.2. Tully-Fisher Distance to Dwl 

In order to determine the distance to Dwl, we chose to apply the IRTF relation 
because of its lower intrinsic dispersion and the reduced extinction in the if-band, with 
the additional advantage that the extinction and the distance are estimated at the same 
wavelength. The IRTF relation was pioneered by Aaronson et al. (1979). However, we 
chose to employ the IRTF relation calibrated by Freedman (1990) using local galaxies with 
Cepheid based distances, and the relations of Peletier & Willner (1993) calibrated relative 
to the distance of the Ursa Major galaxy cluster. 

Freedman's (1990) calibration for the IRTF relation leads to the following expression: 
H^ 5 = -10.26(±0.49)(log AV 20 (0) - 2.5) - 21.02(±0.08), where AV 20 (0) is the inclination 
corrected 20% level Hi velocity profile width in km s -1 . H^q 5 * s the absolute ff-band 
magnitude within a circular aperture with diameter A, for which log(A/D ) = —0.5, 
with D being the 5-band isophotal diameter at fiB,o = 25 mag arcsec~ 2 . Using 
AV 2 o(0) = 259.0 ± 0.5 km s -1 , which is Burton et al.'s (1996) value corrected for the 
inclination of the galaxy, the above expression predicts an absolute iJ-band magnitude 
if***. = _20.13 ± 0.21 mag. 

Although £>-band surface photometry for Dwl is not available, we can make use of the 
intrinsic integrated colors for Sb-Sc galaxies B — H = 3.28 ±0.14 (de Jong 1996) to estimate 
D . The Hb,o — 25 mag arcsec" 2 isophote corresponds to an ff-band surface brightness 
(corrected for extinction) of /1h,o — 21.72 ± 0.14 mag arcsec -2 . Taking into account the 
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if-band extinction this is equivalent to an observed (not corrected for extinction) value 
of hh = 22.19 ± 0.17 mag arcsec -2 . As can be seen from Figure 2 (upper panel), this 
value exceeds the boundaries of the if -band mosaic. However, the surface brightness 
profile can be easily extrapolated. We fitted an exponential disk oc e~ r ^ Td , where r 
is the semi-major axis and r d is the disk scale length) to the surface brightness profile 
from a radial distance of 2 arcmin outwards where the bulge contribution in negligible, and 
estimated a value of D = 8.5 ±0.8 arcmin. The apparent //-band magnitude (not corrected 
for extinction) for a circular aperture with diameter of A_ .5 = 2.7 ± 0.2 arcmin is then 
m(H^ 5 ) = 8.96 ± 0.12 mag, which provides a distance modulus of (m — M) = 28.62 ± 0.26 
and a distance of d — 5.3lo;g Mpc. 

The IRTF relation was initially calibrated for circular apertures because only 
single-pixel detectors were used at the time (Aaronson, Huchra, & Mould 1979). 
Naturally, one would expect a transition to elliptical apertures to reduce the internal 
dispersion of the IRTF relation because they correct for the galaxy inclination. Peletier 
& Willner (1993) studied the problem in detail and reported no significant change in 
the calibration of the IRTF relation when elliptical apertures were used. A possible 
explanation is that the higher internal absorption, which increases with inclination, may 
cancel out the projection effect. Peletier & Willner (1993) used elliptical apertures and 
obtained the following calibration for spiral galaxies in the Ursa Major galaxy cluster: 
log AV2o(0) = -0.085(ir b S 5 e + 30.95 - 9.0) + 2.603. We used a distance modulus to 
Ursa Major Cluster of (m - M) = 30.95 ± 0.17 mag (Pierce & Tully 1988). This gives 
H_ D S £ = -19.72 ± 0.52 mag. However, we are still left with the problem of extrapolating 
the observed total luminosity profile out to D . The measured total (not corrected for 
extinction) //-band magnitude within an elliptical aperture with major axis 2.7 ±0.2 arcmin 
and axial ratio 1.56 (Loan et al. 1996) is m(H^ 5 ) = 9.19 ± 0.43 mag. Correcting for the 
reddening, we obtained a distance modulus of (m — M) = 28.44 ± 0.69 mag and distance 
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ci = 4.9±i;| Mpc. The intrinsic spread of Peletier & Willner (1993) IRTF relation and 
the uncertain distance to the Ursa Major Cluster account for the increased errors of this 
estimate. 

We can now use our reddening estimate to correct the optical / and i?-band photometry 
of Loan et al. (1996), and use the optical Tully-Fisher relations to obtain another distance 
estimate. The two reddening laws discussed in the preceding section predict the same 
Aij An ratio within a few percent. Adopting Mathis (1990) extinction law, we find I- and 
.B-band extinctions of Aj — 1.3 ± 0.3 mag, and As = 3.6 ± 0.8 mag lower than As = 4.3 mag 
in Loan et al. (1996). This puts Dwl at an average distance of 5. 51^7 Mpc. 

These distance determinations place Dwl behind NGC 1560 at 3.5 ± 0.7 Mpc, 
UGCA 105 at 3.8 ± 0.9 Mpc (Krismer, Tully, & Gioia 1995), and Maffei 1 at 4.2 ± 0.5 Mpc 
(Luppino & Tonry 1993). 

4. Conclusions 

We have obtained deep near-infrared imaging of the highly obscured galaxy Dwl. The 
observed infrared colors were used to obtain a very accurate estimate of the extinction 
in the iJ-band, Ah = 0.47 ±0.11 mag. This value was confirmed by the optical - near 
infrared color excesses, and is close to the estimate based on the IRAS 100 /im flux (Loan 
et al. 1996). Our approach is more reliable than previous works in that we did not make 
any additional assumptions for the relation between A v and H I, or the IRAS 100 /iin 
emission. In addition, the IRTF relation allowed us to estimate both the reddening and 
distance at the same wavelength range. This makes our results largely independent of the 
choice of the reddening law, with the additional advantage that the IRTF shows a smaller 
dispersion than its optical counterpart. Finally, the infrared reddening estimates are more 
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reliable than those in the optical because infrared colors of spiral galaxies show a smaller 
intrinsic dispersion, and are less sensitive to the history of star formation than optical colors 
(Vazdekis et al. 1996). 

The IRTF relation (Freedman 1990) yielded a distance of d = 5.3^Q gMpc which places 
Dwl at the far end of the IC 342/Maffei 1 & 2 group. We also confirmed that Dwl has an 
angular diameter greater than 7 arcmin, larger than the originally measured value of 4.2 
arcmin. 



During the course of this work VDI and AA-H were supported by the National 
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Figure Captions 

Figure 1. — f/-band mosaic of Dwl displayed on a logarithmic scale. The orientation is 
north up, east to the left. The field of view is 4.7' x 4.9'. 

Figure 2. — Upper panel: Observed surface brightness profiles in J, H and K s as a 
function of the semi- major axis. The dashed line is our exponential disk fit (see text). 
Middle panel: Radial distribution of the J — H, H — K s and J — K s colors. The straight 
lines represent a linear fit to the color gradients as a function of the semi-major axis (see 
text). Bottom panel: Observed total magnitude as a function of the semi-major axis in J, 
H and K s . In all three panels, the vertical bars represent 3a errors. 
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Table 1. Log of the observations. 

Detector Filter t cxp Date Conditions seeing 

(s) (arcsec) 



256 x 256 J 600 Nov 5 1998 photometric 1.2 

J 780 Jan 1 1999 non-photometric 1.5 

H 1800 Dec 30 1998 non-photometric 1.1 

K s 720 Nov 5 1998 photometric 1.1 

1024 x 1024 H 760 Jan 6 1999 photometric 1.0 
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Table 2. Surface brightness in JHK S for Dwl. 



Semi-major 


Hj 




fJ>H 


a{fi H ) 




a(\i K ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1.2 


18.35 


0.14 


17.21 


0.15 


16.80 


0.18 


3.6 


18.87 


0.06 


17.88 


0.06 


17.42 


0.07 


6.0 


19.09 


0.07 


18.05 


0.06 


17.63 


0.07 


8.4 


19.29 


0.07 


18.18 


0.06 


17.74 


0.07 


10.8 


19.39 


0.07 


18.36 


0.06 


17.88 


0.07 


13.2 


19.57 


0.06 


18.50 


0.06 


18.07 


0.07 


15.6 


19.72 


0.06 


18.62 


0.07 


18.20 


0.07 


18.0 


19.88 


0.06 


18.74 


0.07 


18.33 


0.07 


20.4 


19.98 


0.06 


18.87 


0.07 


18.47 


0.07 


22.8 


20.08 


0.06 


18.96 


0.07 


18.56 


0.07 


25.2 


20.18 


0.06 


19.07 


0.07 


18.69 


0.07 


27.6 


20.27 


0.06 


19.15 


0.07 


18.78 


0.07 


30.0 


20.36 


0.06 


19.25 


0.07 


18.88 


0.07 


32.4 


20.40 


0.06 


19.29 


0.07 


18.91 


0.07 


34.8 


20.46 


0.06 


19.34 


0.07 


19.00 


0.07 


37.2 


20.50 


0.06 


19.41 


0.07 


19.06 


0.07 


39.6 


20.55 


0.06 


19.49 


0.06 


19.09 


0.07 


42.0 


20.59 


0.06 


19.55 


0.06 


19.15 


0.07 


44.4 


20.64 


0.06 


19.59 


0.06 


19.19 


0.07 


46.8 


20.70 


0.06 


19.62 


0.06 


19.27 


0.07 


49.2 


20.74 


0.06 


19.69 


0.06 


19.33 


0.07 


51.6 


20.76 


0.06 


19.76 


0.06 


19.38 


0.07 
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Table 2 — Continued 






Semi-major 


LI j <t( LI A Liu 


lj( Liu) 


r 


a(u K ) 


(1) 


(2) (3) (A) 
y^J 


(5) 


(6) 


(7) 


54.0 


20.80 0.06 19.79 


0.06 


19.43 


0.07 


56.4 


20.87 0.06 19.83 


0.06 


19.43 


0.07 


58.8 


20.91 0.06 19.89 


0.06 


19.47 


0.07 


61.2 


19.91 


0.06 






70.8 


20.04 


0.06 






80.4 


20.13 


0.06 






90.0 


20.21 


0.06 






99.6 


20.35 


0.06 






109.2 


20.44 


0.07 






118.8 


20.53 


0.07 






128.4 


20.68 


0.07 






138.0 


20.85 


0.07 






147.6 


21.01 


0.08 






157.2 


21.13 


0.08 






166.8 


21.08 


0.21 






176.4 


21.38 


0.10 







Note. Column (1) is the semimajor axis in arcsec. 

Columns (2), (4) and (6) are the observed surface brightnesses 
expressed in mag arcsec -2 in J, H, and K s . Column (3), (5) and 
(7) are the 3a errors (see text). 
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Table 3. Colors and redenning. 





(J-H) 


(J-K s ) 


(H - K s ) 


Colors a 


1.04 


1.40 


0.36 


Photometric errors 


0.09 


0.08 


0.09 


Fit errors 


0.08 


0.08 


0.08 


Total errors 


0.12 


0.11 


0.12 


Intrinsic colors 


0.73 ±0.02 


0.94 ±0.03 


0.21 ±0.02 


Color excesses 


0.31 ±0.12 


0.46 ±0.11 


0.15 ±0.12 


A H (M90) 


0.51 ±0.20 


0.47 ±0.11 


0.39 ±0.31 


A H (RL85) 


0.51 ±0.20 


0.46 ±0.11 


0.42 ±0.33 


A v (M90) 


2.93 ±1.13 


2.70 ±0.63 


2.20 ± 1.76 


A v (RL85) 


2.90 ±1.12 


2.76 ±0.64 


2.38 ±1.90 


A B (M90) 


3.88 ± 1.50 


3.57 ±0.83 


2.91 ±2.33 


A B (RL85) 


3.84 ± 1.48 


3.65 ±0.85 


3.15 ±2.52 



a Colors (not corrected for extinction) at 1 arcsec radial 
distance as obtained from the fits (see text). 
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2 1 -1 -2 



R.A. offset (arcmin) 

Fig. 1. — if-band mosaic of Dwl displayed on a logarithmic scale. The orientation is north 
up, east to the left. The field of view is 4.7' x 4.9'. 



-20- 




Fig. 2. — Upper panel: Observed surface brightness profiles in J, H and K s as a function 
of the semi-major axis. The dashed line is our exponential disk fit (see text). Middle panel- 
Radial distribution of the J — H, H — K s and J — K s colors. The straight lines represent 
a linear fit to the color gradients as a function of the semi- major axis (see text). Bottom 
panel: Observed total magnitude as a function of the semi-major axis in J, H and K s . In 
all three panels, the vertical bars represent 3a errors. 



